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Abstract. Classes of relativistic symmetries accommodating supersymmetric patterns are consid- 
ered for the Dirac Hamiltonian with axially-deformed scalar and vector potentials. 
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The Dirac equation plays a key role in microscopic descriptions of many-fermion 
systems, employing covariant density functional theory and the relativistic mean-field 
approach. In application to nuclei and hadrons, the required Dirac mean-field Hamilto- 
nian involves a mixture of Lorentz vector and scalar potentials. Recently, symmetries 
of Dirac Hamiltonians with such Lorentz structure were shown to be relevant for ex- 
plaining observed degeneracies in the spectra of nuclei (pseudospin doublets dll]) and 
mesons (spin doublets [2]). Corresponding supersymmetric patterns were identified for 
spherically- symmetric potentials [3]. In the present contribution, we extend these studies 
to Dirac Hamiltonians with axially-deformed scalar and vector potentials [4]. 

The essential ingredients of supersymmetric quantum mechanics [5] are the su- 
persymmetric Hamiltonian M' = {^^^ ^ ) ^ {^Q^ uJ^ ^'^^ charges Q = (^lo)' 
2+ = {Q-)\ which generate the supersymmetric algebra 

[^,e±] = {e±,e±} = o , {q-m+} = ^. a) 

//+ and //_ satisfy an intertwining relation, LH^ = H L, which ensures that if 0+ is an 
eigenstate of //+, then also (^^ = L<^+ is an eigenstate of //_ with the same energy. Such 



L* 



H. H. 
good SUSY broken SUSY 



FIGURE 1. Typical spectra of good and broken SUSY. The operators L and connect degenerate 
states in the //+ and H sectors. 



TABLE 1. Classes of symmetries and related supersymmetries of the Dirac Hamiltonian 
with axially-deformed scalar, Vs{p,z), and vector, VV(p,z), potentials |4]. 



Class 


Symmetry 


Supersymmetry 


(I) 


Vs{p,z)+Vv{p,z) = Ao (pseudospin) 


(good SUSY) 


(11) 


Vs{p,z)-Vv{p,z) = So (spin) 


(broken SUSY) 


(III) 


Vs^Vs(z),Vv^Vv{p) 


Vs = Vsiz), Vv = f (good SUSY) 


(IV) 


Vs = Vsip),Vv = Vv{z) 


Vs = f, Vv= Vv{z) (good SUSY) 



pairwise degeneracy persists, unless L(^+ vanishes or produces an unphysical state, {e.g. , 
non-normalizable), as occurs for the ground state when the supersymmetry (SUSY) is 
exact. The SUSY partner Hamiltonians and H are therefore isospectral in the sense 
that their spectra consist of pairwise degenerate levels and a possible non-degenerate 
single state in one sector, when the supersymmetry is exact. If all states are pairwise 
degenerate, the supersymmetry is said to be broken. Typical spectra for good and broken 
SUSY are shown in Fig. 1 . 

The Dirac Hamiltonian, H, for a fermion of mass M moving in external scalar V5 and 
vector Vy potentials is given by H = a ■ p + j8 (M + V5) + Vy, where a,;, j8 are the usual 
Dirac matrices, and we have set ^ = c = 1. When the potentials are axially-symmetric, 
i.e., independent of the azimuthal angle 0, Vsy = V5,y(p,z), p = a/jc^ +y-^, then the 
z-component of the angular momentum operator, J^, commutes with the Hamiltonian, 
and its eigenvalues Q. are used to label the Dirac wave functions 



^n(p,(^),z) 



/ ^+(p,z)e'(^-l/2)'^ \ 

^-(p,z)e'(^+V2)<l' 

//+(p,z)e'(^-l/2)^ 



(2) 



Here g^{p,z) and f^{p,z) are radial wave functions for the upper and lower compo- 
nents respectively. The Dirac equation, //^ = E^', leads to four coupled partial differen- 
tial equations for these functions. Their solutions are greatly simplified in the presence 
of symmetries. In what follows we discuss four classes of relativistic symmetries, listed 
in Table 1, and possible supersymmetries within each class 

The symmetry of class I, referred to as pseudospin symmetry, occurs when the sum 
of the scalar and vector potentials is a constant, Vs{p,z) +Vv{p,z) — Aq. The symmetry 

generators, 5,, commute with the Dirac Hamiltonian and span an SU(2) algebra [0,0] 

Si) Up= . (3) 

Here f,- = a,/2 are the usual spin operators, defined in terms of Pauli matrices. The Dirac 
eigenfunctions in the pseudospin limit satisfy 

5,^[f) = A^?^ A = ±1/2 (4) 
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FIGURE 2. Grouping of deformed shell-model states [A^ = 4,«3,A]n, exhibiting a pattern of good 
SUSY relevant to the pseudospin symmetry limit. and «3 are harmonic oscillator quantum numbers. A 
is the pseudo-orbital angular momentum projection along the symmetry z-axis. 



and form degenerate SU (2) doublets. Their wave functions have the form 



(-1/2) 
fli=A-l/2 







1/2) 
fl2=A+l/2 



ge' 



iA<l> 



i{A+l)(j) 



g e 




(5) 



where A = f2 — /I is the eigenvalue of — S^- The relativistic pseudospin symmetry has 
been tested in numerous realistic mean field calculations of nuclei and were found to 
be obeyed to a good approximation, especially for doublets near the Fermi surface |@, 
[Toll . The dominant upper components of the states in Eq. ([5]), involving g^ and g~, 
correspond to non-relativistic pseudospin doublets with asymptotic (Nilsson) quantum 
numbers [A^, ^3 , A] f2 = A + 1 /2 and [A^, «3 , A + 2] t2 = A + 3/2, respectively. The doublet 
is expressed in terms of the pseudo-orbital angular momentum projection, A = A + 1, 
which is added to the pseudospin projection, /i = ±1/2, to form doublet states with 
^2 = A± 1/2. Such doublets play a crucial role in explaining features of deformed 
nuclei, including superdeformation and identical bands |lO(l . The doublet structure can 

be deduced from the fact that the present Dirac Hamiltonian commutes with and L, 



2{M + Ao-H)S, 



(6) 



while L and anticommute. The operator L connects the two states in Eq. (|5]) and 
= L^L = (M + Aq — H)^. Consequently, near the pseudospin limit, the spectrum, for 
each A 7^ 0, consists of twin towers of pairwise degenerate pseudospin doublet states, 
with Qi = A — 1/2 and Q.2 = A + 1/2, and an additional non-degenerate nodeless state 
at the bottom of the = A — 1/2 tower, whose pseudospin partner state is not a bound 
Dirac eigenstate (^]. Altogether, the ensemble of Dirac states with — ^1 = 1 exhibits 
a supersymmetric pattern of good SUSY, as illustrated schematically in Fig. 2. 

The symmetry of class II, referred to as spin symmetry, occurs when the difference of 
the scalar and vector potentials is a constant, Vs{p,z) — Vv(p,z) = Sq. The discussion 
of spin symmetry is similar to that of pseudospin symmetry with the role of upper and 
lower components interchanged. Specifically, the symmetry group is again SU (2) and 
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FIGURE 3. Grouping of deformed shell-model states [A^ = 4,n3,A]i2, exhibiting a pattern of broken 
SUSY relevant to the spin symmetry limit. A is the orbital angular momentum projection. 



its generators 



Si 
UpS.Uj, 



i=x,y,z 



(7) 



commute with the Dirac Hamiltonian. The Dirac eigenfunctions in the spin limit satisfy 

S.^SJ) = M^o^ M = ±l/2 (8) 
and form degenerate SU (2) doublets. Their wave functions are of the form [[loll 



(1/2) 

£^l=A+l/2 







-1/2) 
£^2= A- 1/2 



/ \ 



(9) 



where A = f2 — /i >Ois the eigenvalue of — S^. The upper components of the two 
states in Eq. ^ form the usual non-relativistic spin doublet with a common radial 
wave function g, an orbital angular momentum projection, A, and two spin orientations 
^2 = A± 1/2. The above relativistic spin symmetry is relevant to the structure of heavy- 
light quark mesons f^. Its characteristic doublet structure can be deduced from the fact 
that the present Dirac Hamiltonian commutes with the operators 5, Q and L 



2(M + Zo+H)S, 



(10) 



while L and anticommute. The operator L connects the two states in Eq. Q and = 
L^L = {M + Eo + H)^. These relations imply that near the spin limit, the spectrum, for 
each A 7^ consists of twin towers of pairwise degenerate spin-doublet states with Qi = 
A — 1/2 and Q2 = A + 1/2. None of these towers have a single non-degenerate state. 
Altogether, the ensemble of Dirac states with Q.2 — ^i = —l exhibits a supersymmetric 
pattern which resembles that of a broken SUSY, as illustrated in Fig. 3. 

The Dirac Hamiltonian has additional symmetries when the scalar and vector poten- 
tials depend on different variables. The symmetry of class III occurs when the potentials 
are of the form V5 = Vs{z) and Vy =Vv{p). In this case, the Dirac Hamiltonian com- 
mutes with the following Hermitian operator 



(11) 



where E,- = g' ^. j ■ The Dirac eigenfunctions satisfy 

R^^fi^) = e^(£) . (12) 
A separation of variables is possible by choosing the Dirac wave function in the form 



xj/(e) _ _±_ 
VP 



( Mi(p)vi(z)e'(^-i/2)<i. ^^ 

M2(P)V2(Z)^^'("+1/2)<|. 

/Mi(p)v2(z)e'(^-i/2)<?. 
V -/M2(p)vi(z)e'(^+i/2)'? / 



(13) 



The Dirac equation then reduces to a set of two coupled first-order ordinary differen- 
tial equations in the variable p, and a separate set in the variable z- The wave functions 
depend parametrically on the quantum number e (fT2)) . which can be continuous or dis- 
crete, and |e| plays the role of a mass for the transverse motion. A particular selection 
of potentials within the symmetry class III, is relevant to the study of electron channel- 
ing in crystals Q/s[z) = 0) [11]. A supersymmetric pattern within the present class is 
obtained for the choice VV(p) = OCv / P and Vs{z) arbitrary. The energy eigenvalues are 

^lp!n = |e|/V'l + «v/K + r)^ (np = 0, 1,2, . . .), with r=\/^^- CCy- For Hp > 1, the 

States ^ are degenerate. For rip = only one state is an acceptable solution, which 
has £ > (assuming ay < 0). For each Q. and £ the spectrum resembles a supersym- 
metric pattern of good SUSY, with the towers (H ) of Fig. 1 corresponding to states 
with £ > (£ < 0). The explicit solvability and observed degeneracies in the present 
case can be attributed to the existence of an additional Hermitian operator 

L = pt3S^iJ,r5[H-t3R,]-^{tp)R,^ (14) 

which commutes with H but anticommutes with R^ (fTT)) . The operator L connects the 
doublet states with {np > 1, ±£) and annihilates the single state with {np = 0,£ > 0). In 
addition, = L'^L = jJ{H^ -Rj) + a^Rl 

The symmetry of class IV occurs when the potentials are of the form Vs = Vs{p) and 
Vv = Vv{z). Its discussion is similar to that of class III, with the p- and z-dependence of 
the potentials interchanged. In this case, the following Hermitian operator 

Rp = [M + Vs{p)]t3-ih5{txp)3 (15) 

commutes with the Dirac Hamiltonian and the Dirac eigenfunctions satisfy 

Rp^^i~^) = gvple) . (16) 



Again, a separation of variables is possible with the choice of wave function, 

in a form similar to that of Eq. (fT3l) . The quantum number £ plays now the role of 

an energy for the transverse motion. A particular selection of potentials within the 



symmetry class IV was encountered in the study of the Schwinger mechanism for 
particle-production in a strong confined field (VV(^) = OCyz) [fl2|] . qq pair-creation in 
a flux tube iVs{p = 0) [[l3l] . and the canonical quantization in cylindrical geometry of a 
free Dirac field iVs{p) = Vv{z) = 0) [ffl . A supersymmetric pattern within the present 
class is obtained for the choice Vs{p) = OCs / p {as <0) and Vy (z) arbitrary. The allowed 

values are e = ±My''l - a|/(np + 7)2 {np = 0, 1, 2. . . .), where 7 = ■sj^'^ + a^. For 

rip > I the eigensolutions with ±e are degenerate, E^^^'^ = E^^^ ^. For rip =0 only one 

state, with e > 0, is an acceptable solution. For each Q. and £ the resulting spectrum 
resembles a supersymmetric pattern of good SUSY. As before, the indicated structure is 
due to an additional Hermitian operator 

L = E3|//,75 [M-E3i?p]--^(E-p)j8i?p| (17) 

which commutes with H but anticommutes with Rp (fT5l) . The operator L connects the 
doublet states with {rip > 1, ±e) and annihilates the single state with {rip = 0,e > 0). In 
addition, = L'^L = jJiRj -M^) + a|i?2 . 

The supersymmetric patterns obtained within each of the four classes of symmetries 
share a common origin, namely, the existence of two conserved and anticommuting 

operators R and L. Here R = 5^, 5^,i?^,i?p, given in Eqs. (l3l).(r7l).(fTTI).(fT5l) and L is 
given in Eqs. (|6l).(fT0l). (fT4l) .(fT7l). for classes (I)-(IV) respectively. The operator R has 
eigenvalues which come in pairs of opposite signs. Denoting the eigenvalue of R by r, 
we can define an operator = R/\r\ with eigenvalues ±1, which provide a "parity" 
label to eigenstates of the corresponding Dirac Hamiltonian, H. The operators L, 
and H satisfy [H , = [H ^ L] = { S^y , L} = 0. The operator L connects eigenstates 
of opposite "parity" and its square, l? = TJL = f{H), is a quadratic function of H. 
Supersymmetric charges {Q±) and Hamiltonian (J^), which satisfy the supersymmetric 
algebra ([T]), can now be constructed as 2± = iL(l t ^r) and = = f{H). The 
operator serves as a Z2-grading operator accompanying the supersymmetric algebra. 
This work is supported by the Israel Science Foundation. 
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